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I.  INTRODUCTION 


The  general  objective  of  the  test  series  presented  in  this  report  was 
to  determine  a  concept  of  a  rarefaction  wave  eliminator  suitable  for  use  on 
the  Ballistic  Research  Laboratory's  (BRL)  2.44  m  shock  tube/blast  simulator. 
The  necessary  parameters  for  the  rarefaction  wave  eliminator  were  to  be 
determined  from  both  scaled  (l/48th)  shock  tube  experiments  and  from  results 
of  predictions  obtained  by  the  NASA-  Ames  one-dimensional  computer  hydrocode.1 

The  study  included  three  general  types  of  rarefaction  wave  eliminators 
(RWEs) .  The  first,  Type  A,  consisted  of  a  solid  reflecting  plate*1  placed 
outside  the  open  end  of  the  test  section  of  the  scaled  shock  tube.  The 
standoff  distance  of  the  RWE  was  ried  as  a  function  of  the  input  shock 
overpressure.  The  second  RWE,  Type  B,  consisted  of  vented  plates3* 4  also 
combined  with  standoff  distances.  The  third  RWE  type  tested.  Type  C,  con¬ 
sisted  of  lattice  strips  across  the  end  of  the  test  section.  No  standoff 
distance  was  used  for  this  type  of  RWE. 

The  hydrocode  was  utilized  to  simulate  a  RWE  on  a  computer  generated 
1/48  scale  model  of  the  BRL  2.44  m  shock  tube.  This  computational  RWE 
represents  a  solid  plate  with  one  circular  vented  area  in  the  center  of  the 
plate  and  no  standoff.  It  most  closely  resembles  experimental  Type  B, 
without  standoff. 


II.  EXPERIMENT 

This  section  describes  the  shock  tube  used  for  the  tests,  the  associated 
instrumentation,  and  a  description  of  the  three  general  types  of  rarefaction 
wave  eliminator  models  tested. 

A.  Modified  5.08  cm  Inside  Diameter  S^ock  Tube 

The  BRL  standard  calibration  shock  tube5  was  modified  in  both  driver 
and  test  section  lengths  so  as  to  model  at  l/48th  scale  the  BRL  2.44  m 


j 

Andrew  Mark,  "Computational  Design  of  Large  Scale  Blast  Simulators,"  AIAA 
19th  Aerospace  Sciences  Meetinq ,  January  12  -  15,  1981,  St.  Louis,  MO. 

2 

J.  J.  Yagla,  and  others,  "The  Waves  System  in  Explosively  Driven  Conical 
Shock  Tubes,  "  Second  MABS,  DNA  22? 5p,  The  Naval  Weapons  Laboratory, 
Dahlgren,  VA,  2-5  November  1970. 

3 

J.  R.  Crosnier  and  J.  B.  G.  Monsac,  "Large  Diameter  High  Performance  Blast 
Simulators, "  Fifth  MABS,  Stockholm,  May  25  -  26,  1977. 

4 

W.  Raverdings,  "Preliminary  Design  of  the  2.0  m  Diameter  FOE  -  Driven  Blast 
Simulator,"  Seventh  MBS,  Medicine  Hat,  Alberta,  Canada,  15-17  July  1981. 

George  A.  Coulter,  "Dynamic  Calibration  of  Pressure  Transducers  at  the  BRL 
Shock  Tube  Facility,"  BRL  Memorandum  Report  1S45,  May  1967  (AD  S54508). 
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shock  tube.6  See  Figure  1  for  a  sketch  of  this  shock  tube.  Figure  2  shows 
a  sketch  of  the  5.08  cm  shock  tube  as  modified  to  give  the  l/48th  scale 
desired. 

The  shock  tube  was  operated  in  an  air-air  mode  with  the  downstream  test 
section  end  modified  for  the  RWEs.  The  different  types  of  the  rarefaction 
wave  eliminators  were  tested  for  effectiveness  as  compared  with  the  open  end 
test  results.  The  diaphragms  were  chosen  to  burst  at  the  driver  pressure 
necessary  to  cover  the  desired  input  shock  overpressure  range  of  25  (3.6  psi) 
to  125  (18  psi)  kPa.  This  range  would  include  the  most  useful  range  of  tests 
at  the  BRL  2.44  m  shock  tube. 

Four  transducer  stations  were  located  at  the  distances  from  the  diaphragm 
section  as  shown  on  Figure  2.  These  stations'  locations,  when  scaled  up, 
corresponded  to  the  present  test  stations  for  the  full-sized  shock  tube 
(2.44  m  diameter).  The  pressure  transducers  and  associated  electronics  are 
described  in  the  next  section. 

B.  Instrumentation 


A  schematic  of  the  data  acquisition-reduction  system  is  shown  in 
Figure  3.  Piezoelectric  quartz  transducers  were  used  in  the  shock  tube  test 
section  to  monitor  the  shock  wave  interactions.  The  effects  on  the  wave  sys¬ 
tem  caused  by  the  different  rarefaction  wave  eliminators  were  evaluated  for 
effectiveness  from  the  pressure  transducer  results. 

The  transducers  were  coupled  through  a  power  supply  and  data  amplifiers 
to  a  digitizing  oscilloscope.  Ons? t?  comparisons  of  the  results  were  made 
directly  from  the  oscilloscope.  Final  data  processing  was  completed  with  the 
computer,  printer,  and  plotter.  Tables  and  plots  of  pressure- time  records  for 
the  various  test  stations  are  presented  for  comparison  of  the  results. 

C.  Types  of  Rarefaction  Wave  Eliminators 

Three  general  types  of  rarefaction  wave  eliminators  (RWE)  were  modeled 
for  the  present  test  series.  They  are  shown  in  Figure  4  and  listed  in  the 
data  tables  as  Types  A,  B,  and  C. 

The  least  compli-uted  RWE,  Type  A,  consisted  simply  of  a  solid  circular 
reflecting  plate  spaced  al  varying  distances  from  the  open  end  of  the  shock 
tube  test  section.  The  second  type  of  RWE,  Type  B,  consisted  also  of  a  cir¬ 
cular  plate  with  a  standoff  distance,  but  the  plate  was,  in  addition,  vented. 
Variations  of  the  shape  of  the  vented  opening  were  tested.  A  single  rectan¬ 
gular  hrle,  a  single  circular  hole,  and  a  number  of  circular  holes  were  tried 
for  a  given  standoff  distance.  Both  shape  and  percentage  RWE  opening  were 
varied  to  determine’the  effectiveness  of  the  RWEs.  The  third  type  of  RWE,  Type  C, 


ft 

Brian  P.  Bertrand ,  "BRL  Dual  Shook  Tube  Facility, "  BRL  Memorandum  Report 
2001 ,  August  1969  (AD  693264). 
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1.  DIAMETER  NOT  TO  SCALE 

2.  DIMENSIONS  IN  METRES 


Figure  1.  Sketch  of  BRL  2.44  m  shock  tube 
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Figure  2.  Sketch  of  modified  5.08  cm  inside  diameter  shock  tube 
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Figure  3.  Schematic  of  data  acquisition-reduction  system 
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Figure  4.  Types  of  rarefaction  wave  eliminators 


16 


SOLID 

REFLECTING 

C 

W/STANDOFF 


VENTED 

PLATE 

W/STAND-OFF 


LATTICE, 
NO  STAND-OFF 


was  a  very  simple  version  of  one  used  at  the  French7  blast  simulator  facility 
at  Gramat,  France.  The  programable  moving  Venetian  blind  type  lattice  used 
there,  although  effective,  was  thought  not  to  be  cost  effective  for  use  at  the 
BRL  2.44  m  shock  tube.  Instead  a  simple  fixed  lattice  RWE  was  used  as  a  sub¬ 
stitute  in  the  present  tests.  Lattice  slats  were  removed  to  vary  the  venting 
area. 


III.  RESULTS 

The  results  are  summarized  in  the  data  tables  and  with  pressure- time 
records  to  illustrate  the  effectiveness  of  the  three  general  types  of  rare¬ 
faction  wave  eliminators  tested. 

A.  Data  Table 


The  pertinent  shot  parameters  are  listed  in  Table  1  for  the  types  of 
RWEs  tested  and  the  input  shock  pressure  level  used.  Table  1  also  summarizes 
the  parameters  pertinent  to  the  individual  RWEs,  such  as  standoff  distance 
and  ratio  of  vented  opening  to  total  shock  tube  cross  section. 

Presented  also  are  the  standoff  distances  necessary  to  eliminate  the 
rarefactions  as  completely  as  possible  for  a  Type  A  RWE  for  the  range  of 
input  pressures  from  23  to  124  kPa.  The  vented  area  listed  in  all  cases. is 
the  effective  area.  Obstruction  of  fasteners  used  to  mount  the  RWE  was  . 
subtracted  out  of  the  vent  area.  The  listed  effective  vented  ratio  varied 
from  about  21  percent  to  81  percent  over  the  pressure  range  tested. 

Also  shown  is  a  comparison  of  the  Type  B  vented  plate  RWEs  with  stand¬ 
off  distances  which  were  used  to  furnish  part  of  the  vented  area.  Again 
the  input  pressure  versus  effective  vented  area  appeared  similar  to  the 
solid  plate.  Vent  area  ratios  range  from  43  to  84  percent. 

The  lattice  RWE  was  used  without  a  standoff  distance.  The  braces  used 
to  combine  shots  in  Table  1  indicate  an  interpolation  between  shots  to  find 
the  listed  input  pressure  for  the  given  vented  area  ratio.  The  lattice 
RWE  is  more  cumbersome  and  less  accurate  to  use  for  this  reason.  If  the 
vented  area  is  too  large  for  a  specific  shock  overpressure  the  reflected 
wave  moving  back  up  the  tube  will  be  too  small  and  rarefaction  effects  will 
not  be  completly  eliminated.  Conversely,  a  vented  area  too  small  will 
reflect  a  shock  which  will  be  too  large  and  again  the  loading  on  the  target 
will  not  be  realistic. 

B.  Pressure-Time  Records  without  RWE 


The  first  series  of  records  shown  in  Figures  5-8  illustrate  the 
pressure-time  waveforms  obtained  for  four  side-on  overpressure  levels  at  the 


7 

C.  Languin,  "Simulation  of  Shock  and  Blast  Expansion  Wave  Compensator 
for  2.40  m  diameter  Shock  Tube ,  Description ,  Report  of  Its  Behavior  in 
Response  to  Blast  Waves,"  Centre  d' Etudes  de  Gramat,  Technical  Note 
T-76-29,  May  1976. 
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TABLE  1.  SHOT  PARAMETERS 


Shot  and 
No. 

Station 

No. 

Shock 

Over¬ 

pressure 

kPa 

Type 

RWE 

Standoff 

Distance 

cm 

Vent 

Area 

Ratio 

Comment 

2-82-8 

1 

23.9 

none 

NA 

1.00 

Open  tube 

2 

23.9 

2-82-2 

3 

25 . 2 

4 

23.2 

2-82-7 

1 

54.5 

none 

NA 

1.00 

Open  tube 

2 

54.5 

2-82-1 

3 

55.3 

4 

51.7 

2-82-6 

1 

84.1 

none 

NA 

1.00 

Open  tube 

2 

87.2 

2-82-3 

3 

91.7 

4 

84.1 

2-82-5 

1 

116.7 

none 

NA 

1.00 

Open  tube 

2 

118.6 

2-82-4 

3 

122.0 

4 

110.3 

2-82-24 

2 

25.0 

none 

NA 

1.00 

Open  tube 

3  Stag 

29.0 

2-82-23 

2 

58.6 

none 

NA 

1.00 

Open  tube 

3  Stag 

71.7 

2-82-22 

2 

91.4 

none 

NA 

1.00 

Open  tube 

3  Stag 

128.9 
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TABLE  1.  SHOT  PARAMETERS  (Continued) 


Shock 


Shot  and 

Station 

Over¬ 

Type 

Standoff 

Vent 

No. 

No. 

pressure 

RWE 

Distance 

Area  Comment 

kPa 

cm 

Ratio 

2-82-21 

2 

124.8 

none 

NA 

1 . 00  Open  tube 

3  Stag 

182.7 

2-82-54 

2 

23.4 

none 

NA 

1 . 00  Open  tube 

3 

25.5 

2-82-51 

2 

59.6 

none 

NA 

1 . 00  Open  tube 

3 

62.7 

• 

2-82-52 

2 

91.0 

none 

NA 

1 . 00  Open  tube 

3 

94.5 

2-82-53 

2 

118.5 

none 

NA 

1 . 00  Open  tube 

3 

125.0 

2-82-16 

3 

22.8 

A  -  Solid 

4 

23.4 

Plate 

0.47 

0.212 

2-82-14 

2 

93.1 

A 

3 

95.5 

1.12 

0.698 

2-82-15 

2 

118.9 

A 

3 

124.1 

1.30 

0.811 

2-82-17 

2 

26.3 

A 

3  Stag 

28.3 

0.47 

0.212 

2-82-18 

2 

58.6 

A 

0.90 

0.557 

3  Stag 

72.4 
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TABLE  1.  SHOT  PARAMETERS  (Continued) 


Shot  and  Station 

No .  No  . 

Shock 

Over¬ 

pressure 

kPa 

Type  Standoff 

RWE  Distance 

cm 

Vent 

Area 

Ratio 

Comment 

2-82-19 

2 

90.7 

A 

3  Stag 

122.4 

1.12 

0.698 

2-82-20 

2  Stag 

123.1 

A 

3 

177.9 

1.30 

0.811 

2-82-97 

3 

25.2 

A 

0.67 

0.253 

2-82-100 

3 

39.6 

A 

0.80 

0.415 

2-82-101 

3 

65.1 

A 

1.00 

0.541 

2-81-102 

3 

75.6 

A 

1.10 

0.603 

2-82-104 

3 

100.0 

A 

1.30 

0.727 

2-82-32 

2 

24.5 

B  Hole  Size 

3 

26.5 

2.54  x  3.49  cm 

0.0 

0.438 

2-82-31 

2 

56.9 

B  Hole  Size 

3 

59.6 

2.54  x  3.49  cm 

0.30 

0.540 

2-82-30 

2 

95.8 

B  Hole  Size 

3 

99.3 

2.54  x  3.49  cm 

0.60 

0.729 

2-82-29 

2 

123.4 

B  Hole  Size 

3 

129.6 

2.54  x  3.49  cm 

0.75 

0.820 

20 


n  - - - - -  •  . . . . .  — - — — — _ _ 

i 

I 

i 

TABLE  1.  SHOT  PARAMETERS 

(Continued) 

t 

' ! 

i  ■ 

l 

i 

Shot  and 
No. 

Shock 

Station  Over- 

No.  pressure 

kPa 

T>Te 

RWE 

Standoff 

Distance 

cm 

Vent 

Area  Comment 

Ratio 

J 

2-82-25 

2 

25.0 

B  Hole  Size 

| 

3  Stag 

30.0 

2.54  x  3.49cm 

0.0 

0.438 

j 

2-82-26 

2 

57.9 

B  Hole  Size 

3  Stag 

71.7 

2.54  x  3.49  cm 

0.30 

0.540 

2-82-27 

2 

94.5 

B  Hole  Size 

3  Stag 

127.6 

2.54  x  3.49  cm 

0.60 

0.729 

2-82-28 

2 

125.4 

B  Hole  Size 

3  Stag 

179.3 

2.54  x  3.49  cm 

0.75 

0.820 

- 

2-82-34 

2 

25.3 

B  Hole  Size 

3 

26.5 

3.33  cm  Dia. 

0.0 

0.431 

2-82-35 

2 

58.6 

B  Hole  Size 

3 

60.3 

3.33  cm  Dia. 

0.30 

0.533 

2-82-36 

2 

91.7 

B  Hole  Size 

3 

95.1 

3.33  cm  Dia. 

0.60 

0.722 

2-82-37 

2 

121.3 

B  Hole  Size 

J 

t 

3 

122.4 

3.33  cm  Dia. 

0.75 

0.815 

2-82-49 

2 

25.5 

B  Hole  Size 

* 

3  Stag 

29.6 

3.33  cm 'Dia. 

0.0 

0.431 

\ 

\ 

\ 

2-82-48 

2 

56.5 

B  Hole  Size 

3  Stag 

71.0 

3.33  cm  Dia. 

0.30 

0.533 
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TABLE  1.  SHOT  PARAMETERS  (Continued) 


Shot  and 
No  . 

Station 

No. 

Shock 

Over¬ 

pressure 

kPa 

Type 

RWE 

Standoff 

Distance 

cm 

Vent 

Area 

Ratio 

Comment 

2-82-47 

2 

93.1 

B 

Hole  Size 

3  Stag 

129.6 

3. 

33  cm  Dia. 

0.60 

0.722 

2-82-46 

2 

125. S 

B 

Hole  Size 

3  Stag 

186.2 

3. 

,33  cm  Dia. 

0.75 

0.815 

2-82-41 

2 

24.8 

B 

10  Holes 

3 

26.2 

1 

cm  Dia. 

0.0 

0.452 

2-82-40 

2 

58.6 

B 

10  Holes 

3 

59.3 

1 

cm  Dia. 

0.30 

0.554 

2-82-39 

2 

95.3 

B 

10  Holes 

3 

100.0 

1 

cm  Dia. 

0.60 

0.743 

2-82-38 

2 

120.0 

B 

10  Holes 

3 

123.4 

1 

cm  Dia. 

0.75 

0.836 

2-82-42 

2 

25.2 

B 

10  Holes 

3  Stag 

29.6 

1 

cm  Dia. 

0.0 

0.452 

2-82-43 

2 

54.5 

B 

10  Holes 

3  Stag 

68.3 

1 

cm  Dia. 

0.30 

0.554 

2-82-44 

2 

98.3 

B 

10  Holes 

3  Stag 

138.6 

1 

cm  Dia. 

0.60 

0.743 

2-82-45 

2 

124.1 

B 

10  Holes 

3  Stag 

184.8 

1 

cm  Dia. 

0.75 

0.836 

22 


TABLE  1. 


SHOT  PARAMETERS  (Continued) 


Shot  and  Station 
No .  No. 


2-82-56  2 

3 

2-82-57  2 

3 

2-82-58  2 

3 

2-82-59  2 

3 

2-82-60  2 

3 

2-82-61  2 

3 

2-82-62  2 


2-82-103  3 


2-82-109  3 


Shock 

Over¬ 

pressure 

kPa 

Type 

RWE 

Standoff 

Distance 

cm 

Vent 

Area  Comment 
Ratio 

24.7 

C  Lattice 

26.7 

5  Slats 

0.0 

0.465 

45.0 

C  Lattice 

47.6 

5  Slats 

0.0 

0.465 

58.6 

C  Lattice 

61.7 

2  Slats 

0.0 

0.716 

105.5 

C  Lattice 

111.7 

2  Slats 

0.0 

0.716 

92.4 

C  Lattice 

95.8 

1  Slat 

0.0 

0.868 

103.4 

C  Lattice 

108.9 

1  Slat 

0.0 

0.868 

124.0 

C  Lattice 

130.3 

1  Slat 

0.0 

0.868 

24.7 

C  Lattice 

0.0 

0.465 

5  Slats 

38.7 

C  Lattice 

0.0 

0.465 

5  Slats 
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TABLE  1.  SHOT  PARAMETERS  (Continued) 


Shock 

Shot  and 

Station 

Over¬ 

Typc 

Standoff 

Vent 

No. 

No . 

pressure 

RWE 

Distance 

Area  Comment 

kPa 

Ratio 

2-82-110 

3 

63.6 

C  Lattice 

0.0 

0.597 

4  Slats 

2-82-113 

3 

83.7 

C  Lattice 

0.0 

0.705 

2  Slats 

2-82-116 

3 

171.5 

C  Lattice 

0.0 

0.879 

1  Slat 

NOTES:  1  - 

Cross  Sectional  Area 

of  Shock  Tube 

:  is  20.26 

2 

cm 

2  - 

Ambient 

Temperature  Range  18.2  • 

22.0  °C 

3  -  Ambient  Pressure,  99.1  -  104.1  kPa 
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records  -  test  section  open,  average 
>  24.1  kPa 


XPk 


foui’  test  stations.  These  correspond  to  the  current  test  stations  (multiply 
time  by  48)  in  the  BRL  2.44  m  shock  tube.  The  effect  of  the  two  types  of  rare¬ 
faction  wave  can  be  traced  from  station  to  station.  The  first  rarefaction 
following  the  flat  part  of  the  pressure  records  is  the  one  propagated  upstream 
from  the  open  end  of  the  test  section.  The  second  rarefaction  arrives  from 
the  closed  end  of  the  driven  section  after  about  10  ms.  Both  rarefaction  waves 
degrade  the  positive  wave  dui'ation  and  significantly  change  the  flow  conditions 
behind  the  shock  front. 

Test  Station  3  (87  in  the  2.44  m  shock  tube)  was  of  the  most  interest  since 
the  large  tube’s  test  platform  is  located  at  this  point.  Therefore,  stagnation 
records  were  also  taken  in  the  5.08  cm  shock  tube  corresponding  to  the  test 
platform's  location.  Figure  9  shows  these  records  from  Station  3  over  the 
input  pressure  range  of  the  test. 

Figures  10-13  (bottom  traces)  show  records  of  subtractions  of  the  side-on 
pressure  from  the  stagnation  pressure.  This  difference  is  a  measure  of  the 
dynamic  load  (compressible  Q)  available  to  a  test  target.  Figure  11  (bottom 
trace)  illustrates  the  efiects  of  the  two  rarefaction  waves  on  the  drag  load. 

The  dynamic  pressure  increases  rapidly  at  the  arrival  of  the  open-end  rare¬ 
faction  wave~4  ms  and  decreased  upon  the  arrival  of  the  closed-end  rarefaction 
wave  ^9  ms.  At  this  input  pressure  level  the  drag  load  on  the  face  of  the  test 
target  will  increase  to  over  four  times  the  desired  load.  The  test  target  might,, 
therefore,  be  over-ranged,  making  any  given  test  on  it  invalid. 

C  Pressure-Time  Records  with  RWE 

The  addition  of  a  RWE  on  the  open-end  of  the  shock  tube  provides  ° 
reflected  shock  and  a  rarefaction  wave.  As  the  initial  shock  wave  reaches  the 
RWE,  both  the  reflected  shock  which  comes  from  the  solid  structure  of  the  RWE, 
and  the  rarefaction  wave  which  comes  from  the  opening  in  the  RWE,  begin  to  travel 
up  the  tube  against  the  flow.  By  adjusting  the  opening  of  the  RIVE  the  strength 
of  both  the  reflected  shock  and  the  rarefaction  wave  can  be  tailored  to  minimize 
the  flow  perturbation  in  the  test  section.  Figures  14-17  show  records  taken  from 
Station  3  with  a  Type  A  solid  plate  RWE  in  place  at  the  end  of  the  test  section. 

A  comparison  of  Figure  15  with  Figure  11  illustrates  the  smoothed  drag  loading 
and  more  correct  input  load  for  a  test  of  a  drag  target,  as  well  as  a  longer 
duration  shock  wave. 

Figures  18  -  21  show  the  same  range  of  tests  and  subtractions  of  records 
for  a  Type  B,  vented  plate  RWE.  The  pressure-time  records  do  not  look  sig¬ 
nificantly  different  from  these  ootained  using  the  Type  A,  a  solid  plate  RWE. 

Both  appear  to  be  equally  effective  in  eliminating  the  effects  of  the  rare¬ 
faction  wave. 

The  lattice  RWE,  Type  C,  was  found  to  be  difficult  to  adjust  to  just 
eliminate  the  rarefactions.  Accordingly,  no  stagnation  pressure-time  records 
were  taken  with  the  lattice  RWE.  If  the  reflecting  area  of  the  RWE  is  too 
small  then  the  reflected  wave  will  be  too  small  in  magnitude  and  will  appear 
as  shown  in  Figure  22  on  Shot  58.  Conversely,  if  the  reflecting  area  is  too 
large  then  the  reflected  wave  will  be  too  large  and  will  appear  as  shown  in 
Shot  57  in  Figure  22.  The  venting  area  at  the  end  of  the  shock  tube  must  be 
matched  with  the  overpressure  reaching  the  end  of  the  tube  in  order  to  elim¬ 
inate  the  effects  of  the  rarefaction  wave. 

Written  material  continued  on  page  46. 
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TEST:  RWE,  2  IN  ST 


TEST:  RWE  2  IK  ST 


Figure  10.  Stagnation  pressure  minus  side-on  -  test  section 
open,  average  input  pressure  24,2  kPa 
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TEST:  RVE  2  IN  ST 


TEST:  RVE,  2  IN  ST 


Figure  11.  Stagnation  pressure  minus  side-on  -  test  section 
open,  average  input  pressure  59.0  kPa 
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TEST:  RUE  2  IN  ST 


TEST:  RUE,  2  IN  ST 


Figure  12.  Stagnation  pressure  minus  side-on  -  test  section 
open,  avei'age  input  pressure  91.2  V;Pa 
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TEST:  RUE,  2  IN  ST 


Figure  13.  Stagnation  pressure  minus  side-on  -  test  section 
open,  average  input  pressure  121.1  kPa 
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TEST;  RWE,  2  IN  ST 


TEST:  RWE,  2  IN  ST 


TEST:  RVE,  2  IN  ST 


Figure  14.  Stagnation  pressure  minus  side-on  -  Type  A  solid 
plate  RWE,  average  input  pressure  24,5  kPa 
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Figure  22.  Pres sure- time  records  for  Type  C  lattice  RWE 


TEST:  RUE,  2  IN  ST 
SHIT!  61 
STATION:  2 


18  15  28  -25 

TDE.  NSEC 


TEST •  RUE,  2  IN  ST 
SHOT.  61 
STATEN:  3 


18  15  28  25 

TBt,  NSEC 


TEST:  RUE.  2  IN  ST 
SHOT.  62 


TEST:  RUE.  2  IN  ST 


ressure 


4  5 


-time  records  for  Type  C  lattice  RWE 


A  comparison  of  records  taken  when  various  vented  RWEs  were  in  place 
with  one  taken  when  the  solid  plate  RWE  was  in  place  is  shown  in  Figure  23 
for  an  intermediate  input  shock  pressure.  The  vented  plate  RWE  was  varied; 
making  use  of  a  single  rectangular  opening,  a  single  circular  opening,  and 
ten  smaller  holes.  The  vented  opening  was  equal  for  all  RWEs  used  to  obtain 
the  results  shown  in  Figure  23.  All  RWEs  appear  to  be  equally  effective. 

The  lattice  RWE  was  not  compared  here  because  of  the  previous  difficulties 
in  its  use. 

Additional  pressure-time  records  for  the  .  different  RWEs  are  shown  in 
Appendix  A. 


IV.  NASA-AMES  COMPUTER  CODE  PREDICTIONS 


A.  Introduction 


Utilizing  a  mainframe  computer  it  is  possible  to  simulate  complicated 
fluid  flow  by  employing  a  mathematical  algorithm  based  on  either  the 
Lagrangian  or  Eulerian  fluid  mechanics  model.  Detailed  information  describing 
field  variables  may  be  obtained  at  any  spacial  and  temporal  location  which 
is  essentially  the  computational  equivalent  of  a  fluid  dynamics  experiment. 
Although  mathematically  possible  to  perform  this  task  without  a  computer  it 
is  impractical  because  of  the  large  number  of  calculations  required. 

B.  Theory 


The  NASA-Ames  hydrodynamic  code  employed  in  this  study  is  a  one-dimen¬ 
sional,  adiabatic,  inviscid,  Eulerian  computer  algorithm  which  was  written  by 
Dr.  Andrew  Mark* and  modified  by  Mr.  Klaus  Opalka^f  the  Ballistic  Research 
Laboratory. 

The  Euler  equations  for  conservation  of  mass,  momentum,  and  energy  per 
unit  volume  are  solved  in  differential  form  for  the  field  variables:  density, 
pressure,  temperature,  total  energy,  and  a  one-dimensional  component  of  flow  9 
velocity  using  finite  differencing  formulations  attributed  to  Beam  and  Warming.' 
The  ideal  gas  equation  of  state.  Equation  1,  and  the  Euler  equations.  Equation  2, 


Q 

Patrick  J.  Roache,  Computational  Fluid  Mechanics .  Hermosa  Publishers , 
Albuquerque,  NM,  1972 ,  pp  2  -  12,  204  -  286. 

9 

R.  M.  Beam  and  R.  F.  Harming,  "An  Implicit  Factored  Scheme  for  the 
Compressible  Navier  Stokes  Equations,"  AIAA  Journal ,  Volume  l 6,  No.  4, 

April  1978,  pp  393  -  402. 

*  Private  communication  with  Dr.  Andrew  Mark,  Ballistic  Research  Laboratory, 
Jan  1982. 

**Private  communication  with  Mr.  Klaus  Opalka,  Ballistic  Research  Laboratory, 
Jan  1982. 
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TEST:  RUE,  2  IX  ST 


TB€,  MSEC 


A.  Solid  plate  RWE  . 
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B.  Rectangular  vented  RWE. 


Figure  23.  Comparison  of  the  effectiveness  of  the  Type  A  solid 
plate  RWE  with  the  Type  B  vented  RWE 
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TEST  i  WE.  2  IN  ST 


C.  Circular  vented  Rif!-. 


TEST:  RUE.  2  IN  ST 


D.  Multiple  circular  vented  RWI-. 


Figure  23.  (continued)  Comparison  oC  the  effectiveness  of  the  Type  A 
solid  plate  RWE  with  the  Type  B  vented  RWE 


(1) 


are  applied  to  a  straight  shock  tube  with  a  RVJE  placed  at  the  open  end. 


p  =  (y-1)  (e-1^2), 

+  i_(puA)=  o 
ot  ox  - 


3  (ppA)  +  _3|(pu2  +  p)A  ] 
3t  3x 


and 


(2-a) 

(2-b) 


3(eA)  ,  3[uA(e+p)l 

“TF"  +37  -  °’  (2-c) 


where  p  =  pressure,  y  is  the  ratio  of  specific  heats,  e  =  total  energy/volume,, 
p  =  density,  u=  flow  velocity,  A  =  tube  cross-sectional  area,  t  =  time, 
and  x  =  distance. 


The  initial  conditions  are  normalized.  Variables  at  the  closed  end  of 
the  tube,  i.e.,  the  compression  chamber  backwall,  are  computed  using  image 
points.  Variables  of  the  open  end  are  calculated  by  using  backward 
differencing. 

Independent  variables  (x,t)  are  transformed  into  a  computational  grid 
and  the  governing  equations  are  solved  at  one-dimensional  spacial  grid  points 
(x)  as  a  function  of  time.  Refer  to  Figure  24,  a  sketch  of  the  computational 
shock  tube.  The  distribution  and  total  number  of  grid  points  are  established 
as  computer  input  parameters.  Typically, the  number  of  spacial  grid  points 
varies  from  100  -  1000  with  200  -  400  providing  adequate  results. 

$ - 3.22 - 5 

$ -  1.73 - % -  1.46 - 5|e-0.83— 5 


Diaphragm 


Grid  Polnls 
(100-1000) 


NOT  TO  SCALE 
Measurements  in  metres 


Figure  24.  Sketch  of  computational  shock  tube. 
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The  computational  grid  may  be  equidistant ly  partitioned  along  the  tube 
length  or  clustered  about  a  specified  location  utilizing  the  hyperbolic 
function  incorporated  in  the  code.  Thus, a  proportionally  large  number  of 
grid  points  may  be  positioned  where  cross-sectional  area  changes  occur, 
i.e.,  at  the  RWE. 

The  RWE  positioned  at  the  open  end  of  the  tube  partially  reflects  the 
shock  wave  which  decreases  the  magnitude  and  effect  of  the  open-end  rarefaction 
wave.  The  hydrocode  models  a  RWE  that  is  a  flat  circular  plate  of  the  same 
diameter  as  the  shock  tube.  In  this  study  the  RWE  has  one  circular  hole  of 
time-invariant  cross-sectional  area  that  allows  for  outflow.  (The  code  also 
provides  for  a  RWE  that  varies  as  a  function  of  time;  this  feature  was  not 
utilized  in  the  present  study.)  The  RWE  percent  of  vented  area  ratio 
(cross-sectional  open  area/ total  cross-sectional  area  X  100)  is  varied  with 
each  computational  run  depending  upon  the  driver  pressure. 

There  are  several  computational  difficulties  with  the  RWE  model  that  do 
not  have  experimental  counterparts.  The  RWE  can  never  be  truly  perpendicular 
to  the  end  of  the  shock  tube.  For  if  it  were,  then  the  last  grid  point  would 
be  dual  valued.  Therefore, the  RWE  is  placed  at  a  slight  angle  and  the  change 
in  tube  area  occurs  over  a  small  distance  x  which  is  also  an  input  parameter. 

This  change  in  cross-sectional  tube  area,  describes  a  parabolic  arc, 

i.e.,  tube  radius  r  is  dependent  on  x, and  A  is  proportional  to  r2. 

The  computational  scheme  is  sensitive  to  the  number  of  grid  points 

within  the  length  of  the  RWE.  Since  is  nonlinear,  a  large  number  of 

grid  points  is  required  to  precisely  define  the  continuous  area  change. 
Computational  results  approach  an  asymptotic  value  if  10  grid  points  are 
placed  within  the  length  of  the  RWE. 

The  imposition  of  an  artificial  length  to  the  RWE  changes  slightly  the 
physical  dimensions  of  the  tube  which  alters  the  normalized  locations  of  the 
test  stations.  This  effect  is  minimal;  the  RWE  length  was  established  as 
1%  of  the  tube  length. 

C.  Results 


The  results  apply  equally  to  the  5.08  cm  and  2.44  m  shock  tubes.  Com¬ 
putational  shock  tube  dimensions  are  identical  to  the  5.08  cm  experimental 
tube.  These  dimensions  are  normalized  by  the  computer  code.  Since  the 
5.08  cm  tube  is  a  l/48th  scale  model  of  the  2.44  m  tube,  the  normalized 
dimensions  apply  to  both  shock  tubes.  For  the  2.44  m  tube, time  of  arrival 
and  shock  duration  must  be  increased  by  a  factor  of  48. 

As  indicated  in  the  theory  section,  the  code  is  sensitive  to  the  number 
of  points  within  the  RWE  length.  By  using  as  input  parameters  404  grid 
points,  a  1%  RWE  length,  and  a  hyperbolic  sine  clustering  factor  8=5, 
seven  grid  points  were  placed  within  the  RWE  length.  Reported  results  are 
within  5%  of  results  achieved  with  10  grid  points  within  the  RWE  length. 

(10  grid  points  provide  a  close  approximation  to  the  asymptotic  value.) 
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Seven  points  were  chosen  to  keep  computer  processing  time  and  output  file 
lengths  within  manageable  limits  while  also  maintaining  accurate  results. 

Figure  25  displays  hydrocode  generated  side-on  and  dynamic  pressure 
(q  =  kfU')  for  an  input  pressure  of  26.8  kPa  at  test  stations  1,  2,  3,  and  4 
without  a  RWE.  Manifestations  of  wave  decay  due  to  both  rarefaction  waves 
are  present.  Stations  1,  2,  3,  and  4  correspond  to  experimental  Stations 
84,  86,  87,  and  88,  respectively, for  the  BRL  2.44  m  shock  tube.  These  results 
are  quite  similar  to  experimental  results  at  a  comparable  imput  pressure. 
Refer  to  Figure  5.  Note  the  zero  time  on  the  experimental  curve  is  the  time 
of  arrival  at  the  station,  while  the  zero  time  on  the  calculations  starts 
at  the  break  of  the  diaphragm.  Direct  computational  and  experimental  com¬ 
parisons  for  various  conditions  are  reported  in  Section  V. 

Figure  26  shows  results  for  62.7,  91.7,  and  122.0  kPa  at  Station  3  which 
corresponds  to  the  primary  test  station  in  the  2.44  m  tube.  Again,  no  RWE 
is  utilized  and  the  rarefaction  wave  effects  are  evident.  Refer  to  Figures 
6,  7,  and  8  for  comparison  with  experiment. 

The  results  displayed  in  Figures  27  through  31  show  the  effects  of  a  RWE 
over  a  large  range  of  input  pressures,  including  input  pressures  used  in 
Figure  26  for  comparison.  All  data  are  computed  at  Station  3.  These  results 
are  reported  in  tabular  form  in  Table  2,  which  shows  the  appropriate  vented 
area  ratios  needed  to  extend  the  shock  wave  duration  as  a  function  of  input 
pressure  at  Station  3. 

TABLE  2.  VENTED  AREA  RATIO  VERSUS  INPUT  SHOCK 
OVERPRESSURE  -  COMPUTER  CODE 


Input  Pressure  Vented  Area 
(kPa)  Total  Area 

2.5  .060 

5.0  .100 

9.8  .150 

21.4  .250 

26.6  .288 

34.3  .339 

40.9  .380 

56.1  .470 

61.4  .500 

71.3  .550 

79.9  .590 

97.2  .663 

119.6  .750 

125.4  .770 

142.4  .820 
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STRTION  2 


STATION  3 


STATION  4 


TIME  CS) 


Figure  25.  Computer  side-on  and  dynamic  pressure-time  records 
26.8  kPa,  no  RWE,  Stations  1,  2,  3,  and  4 


Figure  26.  Computer  side-on  and  dynamic  pressure-time  records 
62.7,  91.7,  and  122.0  kPa,  no  RWE,  Station  3 
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Figure  29.  Computer  side-on  and  dynamic  pressure-time  records, 
40.9,  56.1,  and  61.4  kPa,  RWE,  Station  3 


o 


VENTED  AREA  82.0  X 


Figure  31.  Computer  side-on  and  dynamic  pressure-time  records, 
119.6,  125.4,  and  142.4  kPa,  RWE,  Station  3 


A  secondary  wave*  which  first  appears  at  a  50  percent  vented  area  with  input 
shock  overpressure  of  61.4  kPa  is  a  compression  wave  propagating  upstream. 

This  wave  is  caused  by  the  interaction  of  the  contact  surface  with  the  open 
end  of  the  tube.  The  time  of  arrival  and  magnitude  is  a  function  of  vet. ted 
area  and  shock  overpressure.  See  Figures  30  and  31. 

D.  Discussion 


For  low  to  moderate  pressure  levels  the  RWE  extends  the  duration  of  the 
square  shock  wave  and  reduces  the  dynamic  pressure.  Effects  of  the  open 
end  rarefaction  wave  are  nullified.  At  higher  pressure  the  RWE  extends  the 
duration  of  the  square  wave  and  continues  to  significantly  decrease  the 
dynamic  pressure. 

Figure  32  displays  the  RWE  vented  area  ratio  versus  overpressure  based 
on  the  data  in  Table  2.  This  curve  can  be  used  to  predict  the  effects  of  a 
RWE,  without  standoff,  by  varying  the  circular  opening. 

V.  ANALYSIS 

The  analysis  section  will  treat  the  experimental  results  obtained  from 
the  tests  in  the  5.08  cm  shock  tube,  determine  the  parameters  needed  to 
apply  results  to  the  BRL  2.44  m  shock  tube,  and  compare  experimental  results 
with  predictions  from  the  NASA-Ames  one-dimensional  hydrocode. 

A.  Data  from  5.08  cm  Shock  Tube 


Bertrand*  found  from  experiments  at  the  BRL  57.47  cm  shock  tube  that  a 
larger  reflecting  surface  of  a  test  target  at  the  end  of  the  open  test  section 
had  the  effect  of  a  rarefaction  wave  eliminator.  He  found  that  the  standoff 
distance  and  vented  area  needed  to  eliminate  the  rarefaction  were  related  to 
the  input  shock  overpressure.  Accordingly,  the  data  from  the  present  exper¬ 
iments  were  plotted  as  input  shock  overpressure  as  a  function  of  vented  area 
ratio.  Figure  33  shows  the  data  plotted  in  this  manner  for  the  three  general 
types  of  RWE  tested.  The  asterisks  are  points  from  theory  given  in  Reference  10 
by  Weidermann.  As  noted  in  Table  1  the  vented  plates  with  no  standoff  had  a 
vented  area  ratio  ranging  from  43  to  45  percent,  and  it  can  be  seen  in  Figure 
33  that  those  data  points  do  not  fit  the  trend  of  data  from  RWEs  with  stand¬ 
off  distances.  Figure  34  shows  data  from  the  solid  plate  RWE.  Either  of 
these  curves  (Figure  33  or  34)  may  be  used  to  determine  the  parameters  for  a 
suitable  rarefaction  wave  eliminator  on  the  BRL  2.44  m  shock  tube. 

B.  Parameters  for  Rarefaction  Wave  Eliminator  -  2.44  m  Shock  Tube 


Both  t  e  solid  plate  RWE  and  the  vented  Plate  RWE  appeared  suitable  for 
use  at  the  large  shock  tube.  No  attempt  was  made  to  adapt  the  lattice  type 
of  RWE  to  the  large  shock  tube. 


10 

A.  Weidermann,  "Study  of  Grids  in  Shook  Tubes, "  Proceedings  of  First  Shock  Tube 
Symposium ,  SWE-TM-57-2 ,  Hq  AFSWC,  Kirtland  Air  Force  Base ,  NM,  26-27  Feb  1957. 

*Private  communication  from  Mr. Brian  P. Bertrand,  Ballistic  Research  Laboratory, 
July  1981  and  June  1982. 
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Figure  32.  RWF,  vented  area  ratio  versus  input  shock  overpressure 
-  computer  code 
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HOT:  SOLID  PLATE 


The  solid  plate  needs  only  one  parameter  -  the  standoff  distance  -  to 
be  varied  to  eliminate  the  rarefaction  wave.  The  vented  plate  RWE  needs  two 
parameters  -  the  vented  area  of  the  plate  and  the  standoff  distance. 

Since,  at  present,  the  large  shock  tube  has  an  end  flange  test  stand 
assembly  which  has  a  1.219  x  1.691  m  rectangular  opening,  either  a  solid  or 
vented  RWE  would  be  suitable  using  the  available  hardware.  See  Appendix  B 
for  drawing  details  of  the  flange. 

Calculations  of  standoff  distance*  w,  assumed  twenty-one  bolts/ spacers  of 
4.76  cm  dia.  each.  The  area  obstructed  by  the  bolts  in  the  standoff  zone 
was  subtracted  from  the  total  vented  area  (side  vented  area  plus  area  of  RWE 
plate  hole) .  The  effective  vented  area  ratio  was  used  as  found  from  the 
small  shock  tube  results.  Working  backwards  from  the  given  cross-section 
area,  4.669  (50.26  ft-),  of  the  2.44  m  shocks  tube  and  required  vented  area 

ratio,  the  standoff  was  calculated  fot  each  RWE.  If  a  circular  RWE  is 
assumed,  simple  relationships  may  be  obtained. 

Equations  3-5  below  give  the  results  for  a  solid  plate  RWE. 


Av 


A 


t 


e  R, 


(3) 


the  vented  area  ratio  found  in  Tables  3  and  4  as  functions  of  input  shock 
pressure,  where  Av  is  the  effective  vented  area  and  A  is  the  total  cross 
section  of  the  shock  tube. 


Av  =  (*DW  -  ABolt$  )  -  RAt,  (4) 

the  effective  side  vented  area,  where  the  area  of  the  exposed  bolts 
was  subtracted.  Assuming  twenty-one  bolts.  Equation  4  becomes: 

W  =  0.701R,  m.  (5) 

Using  the  same  method  and  assumptions  from  above,  the  relationship  for 
standoff  for  the  end  flange  assembly's  vented  plate  (2.043  m  ),  is  given 
in  Equation  6. 


W  =  0.701R  -  .307,  m.  (6) 

The  limit  of  R  at  W  -  0  is  0.438  corresponding  to  an  input  of  45  kPa. 

For  values  less  then  this,  the  given  opening  in  the  vented  RWE  plate  will 
have  to  be  reduced  according  to  values  listed  in  the  bottom  half  of  Table  4. 

Figure  35  shows  two  sample  curves  plotted  from  Tables  3  and  4  for  the  RWEs 
used. 
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TABLE  3.  STANDOFF  DISTANCE  VERSUS  INPUT 
SHOCK  OVERPRESSURE  -  SOLID  PLATE  RWE 


Input 

kPa 

Overpressure 

psi 

Vented 

Area 

Ratio 

Standoff 

cm 

Distance 

in. 

0 

0 

0 

0 

0 

5 

0.725 

0.055 

3.8 

1.50 

10 

1.450 

0.116 

8.1 

3.19 

20 

2.900 

0.225 

15.8 

6.22 

30 

4.351 

0.320 

22.4 

8.82 

40 

5.802 

0.400 

28.0 

11.02 

50 

7.252 

0.470 

33.0 

12.99 

60 

8.702 

0.535 

37.5 

14.76 

70 

10.15 

0.593 

41.6 

16.38 

80 

11.60 

0.648 

45.4 

17.87 

90 

13.05 

0.693 

48.6 

19.13 

100 

14.50 

0.735 

51.5 

20.28 

110 

15.95 

0.775 

54.4 

21.42 

120 

17.40 

0.810 

56.8 

22 . 36 

130 

18.85 

0.840 

58.9 

23.19 

140 

20.31 

0.865 

60.7 

23.90 

150 

21.76 

0.885 

62.1 

24.45 

160 

23.21 

0.898 

62.9 

24.76 

170 

24.66 

0.910 

63.8 

25.12 

*Standoff  was  calculated  from  the  vented  area  ratio  based  on  shock  tube 
internal  diameter.  Exposed  bolts/spacers  were  subtracted  from  the  vented 
area.  Twenty-one  bolt/ spacers  were  assumed  -  4.76  cm  4ia.  each. 
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TABLE  4.  STANDOFF  DISTANCE  VERSUS  INPUT 
SHOCK  OVERPRESSURE  -  VENTED  PLATE  RWE 


Input 

Overpressure 

Vented 

Area 

Standoff 

Distance* 

Area 
Mole  in 

of 

Plate 

kPa 

psi 

Ratio 

cm 

in 

m2 

f2 

45 

6.527 

0.437 

0 

0 

2.043 

22.0 

50 

7.252 

0.470 

2.3 

0.91 

60 

8.702 

0.535 

6.8 

2.68 

70 

10.15 

0.593 

10.9 

4.29 

80 

11.60 

0.648 

14.7 

5.79 

90 

13.05 

0.693 

17.9 

7.05 

100 

14.50 

0.735 

20.9 

8.23 

110 

15.95 

0.775 

23.7 

9.33 

120 

17.40 

0.810 

26.1 

10.28 

150 

18.85 

0.840 

28.2 

11.10 

140 

20.31 

0.865 

30.0 

11.81 

150 

21.76 

0.885 

31.4 

12.36 

160 

23.21 

0.898 

32.3 

12.72 

170 

24.66 

0.910 

33.1 

13.03 

0 

0 

0 

0 

0 

0 

0 

5 

0.725 

0.055 

0 

0 

0.257 

2.764 

10 

1.450 

0.116 

0 

0 

0.542 

5.831 

20 

2.900 

0.225 

0 

0 

1.051 

11.31 

30 

4.351 

0.320 

0 

0 

1.494 

16.09 

40 

5.802 

0.400 

0 

0 

1.868 

20.11 

*Standof f 

was  calculated 

from  the 

vented  area 

ratio  for  a 

vented  plate  with 

one  hole  -  1.219  x  1.676  m.  The  ratio  was  based  on  the  internal  diameter  of 
the  shock  tube.  Exposed  bolts/spacers  were  subtracted  from  total  vented  area. 
Twenty-one  bolt/spacers  were  assumed  -  4.76  cm  dia.  each. 


Figure  35.  Predicted  standoff  distance  for  RKEs  on  2.44 
BRL  shock  tube 


C .  NASA-Ames  Computer  Code  Comparisons 


Figures  36  -  41  on  the  following  pages  display  comparisons  between 
experimental  and  computational  results  for  side-on  and  dynamic  pressure 
at  three  pressure  levels  with  and  without  a  RWE.  Results  ar-  for  Station 
3.  Because  the  input  pressure  comparison  cases  are  not  identical,  slight 
differences  in  pressure  levels  arc  evident. 

The  reader  is  advised  to  ohserve  the  general  waveforms.  It  is  apparent 
that  both  experimental  and  computational  results  may  be  used  to  model  a  RWE 
for  the  2.44  m  shock  tube.  There  is  a  difference  between  the  experimental 
and  computational  RWE  vented  areas  used  in  Figures  37,  39,  and  41.  Compar¬ 
isons  of  experimental  (solid  plate)  and  computational  RWE  vented  area  ratios 
as  a  function  of  input  pressure  are  presented  in  Figure  42.  The  two  curves 
are  markedly  different  and  this  distinctness,  which  clearly  is  to  be  expected, 
is  due  to  the  geometrically  different  methods  used  to  obtain  vented  areas. 
Experimentally  a  standoff  distance  is  the  RWF.  mechanism;  this  allows  for 
outflow  normal  to  the  side  of  the  shock  tube.  Computationally,  the  one¬ 
dimensional  hydrocode,  RWE  model  allows  for  outflow  normal  to  the  end  of 
the  tube. 

Although  the  results  are  distinct,  both  methods  of  modeling  a  RWE  are 
useful  and  if  conjoined  provide  an  efficient,  cost-effective  blast  modeling 
technique.  The  experimental  predictions  may  be  used  straightforwardly  to 
forecast  the  appropriate  RIVE  parameters  for  a  specific  shot  in  the  2.44  m 
cube.  Or  to  save  shot  time  and  money  a  first  approximation  to  experimental 
results  may  be  obtained  by  running  the  hydrocode.  Subsequently,  Figure  42 
provides  a  method  for  conversion  from  computational  vented  area  to  exper¬ 
imental  vented  area  applicable  to  the  2.44  m  shock  tube. 

D .  Comparison  of  1/48  Scale  and  Computer  Results  with  Full-Scale  Results 

A  calibration  test  was  conducted  in  the  full-size  tube  with  the 
rarefaction  wave  elimination  standoff  distance  established  from  the  curve 
presented  in  Figure  35  for  a  side-on  overpressure  of  60  kPa.  The  results  of 
this  test  are  shown  in  Figure  43.  The  recorded  side-on  overpressure  was 
64.0  kPa  which  gave  a  higher  reflected  pressure  on  the  RWE  and  therefore  a 
higher  reflected  pressure  moving  back  up  the  tube  as  shown  on  the  records  in 
Figure  43. 

A  computer  run  was  made  to  predict  the  overpressure  versus  time  and  the 
dynamic  pressure  versus  time  at  stations  84,  86,  87,  and  88  in  the  full-size 
(2.44  m)  shock  tube  with  a  rarefaction  wave  eliminator.  These  results  ore 
presented  in  Figure  44.  The  planned  maximum  side-on  pressure  was  64  kPa,  but 
the  rarefaction  wave  eliminator  was  set  for  a  vent  area/ shock  tube  area  ratio 
from  Figure  32  for  a  60  kPa  overpressure  to  match  the  experimental  conditions 
Note  that  the  reflected  wave  from  the  rarefaction  wave  eliminator  is  slightly 
higher  than  the  side-on  pressure  as  shown  in  Figure  43. 
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A  computer  run  was  also  made  for  a  64  kPa  side-on  overpressure  in  the 
full-size  tube  without  a  rarefaction  wave  eliminator.  The  side-on  over¬ 
pressure  and  dynamic  pressure  versus  time  are  presented  for  the  four  stations 
in  Figure  45.  A  comparison  of  Figures  45  and  44  will  show  the  improvement 
in  the  extended  duration  of  the  overpressure  versus  time  and  the  realistic 
dynamic  pressure  when  a  rarefaction  wave  eliminator  is  applied. 

In  Figure  46  a  comparison  is  made  of  the  side-on  overpressure  and 
dynamic  pressure  at  Station  87  for  the  1/48  scaled  shock  tube  model, 
the  computer  code  run,  and  the  full-size  shock  tube,  without  a  rarefaction 
wave  eliminator.  All  pressures  are  normalized  to  an  input  wave  of  62  kPa  for 
direct  comparisons  and  the  time  scale  of  the  record  from  the  model  was 
multiplied  by  48.  The  three  sets  of  records  show  excellent  correlation. 

A  similar  comparison  is  made  in  Figure  47  showing  again  the  correlation 
at  Station  87  between  the  two  methods  used  to  predict  the  effect  of  a 
rarefaction  wave  eliminator  when  placed  to  give  the  proper  vent  area. 


Written  material  continued  on  page  81. 


68 


X  EXPERIMENT 
—  HYDROCODE 


DYNAMIC  PRESSURE 


X  EXPERIMENT 
—  HYDROCODE 


*  K  * 


X  0.015  * 


*  *  '7?'o3^s 


X  ^ 
\  X^  / 


..••••■25--A 
A  A 


A  experiment 

••••  HYDROCODE 


A 

A‘\ 


A  a  A  A 


TIME(S) 


Figure  38.  Comparison  of  experimental  and  computer  side-on  and 
dynamic  records,  62.7  average  kPa,  no  RWE,  Station  3 
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Figure  42.  Comparison  of  experimental  (solid  plate)  and  computational 
predictions  for  rented  area  ratio  versus  input  overpressure 
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STATION  84 


Figure  43.  Overpressure  versus  time  at  four  stations  in  the  2.44 
metre  shock  tube  with  RIVE. 
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Figure  45.  Overpressure  and  dynamic  pressure  versus  time  from  NASA-Ames 
computer  program  for  the  2.44  metre  shock  tube  without  RWE. 
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Figure  46,  Comparison  of  results  from  1/48  scale  model  and  computer  program 
output  with  full-scale  results  at  Station  87  without  RIVE 


79 


PRESSURE (KPRJ  PRESSURE (KPfl) 

-SO.O  0.0  SO.O  100.0  -S0.0  0.0  50.0  100.0  -50.0 


VI.  SUMMARY  AND  CONCLUSIONS 


Three  general  types  of  rarefaction  wave  eliminators  (RWEs)  were 
designed  and  tested  on  a  5.08  cm  I.D.  shock  tube  to  determine  a  suitable 
design  of  RWE  for  the  BRL  2.44  m  shock  tube.  The  small  shock  tube  had  been 
modified  from  a  standard  calibration  shock  tube  to  represent  l/48th  scale 
of  the  larger  shock  tube.  Lengths  of  driver,  test  section,  and  transducer 
locations  were  scaled  to  match  the  larger  shock  tube. 

The  tests  included  parameter  changes  of  input  pressure,  standoff 
distance,  and  plate  opening  shape  and  number  of  openings  in  a  particular 
RWE  plates  used.  Of  the  three  types  tested  -solid  plate,  vented  plate  and 
lattice  strips  -  the  solid  plate  and  vented  plate  types  were  found  to  be 
most  suitable  for  use  on  the  BRL  2.44  m  shock  tube.  It  was  difficult  to  find 
the  exact  vented  area  ratio  needed  when  the  lattice  RWE  was  tried. 

Variations  in  the  shape  of  the  RWE  plate  opening  did  not  seem  to  make 
any  noticeable  difference  in  the  pressure- time  records  obtained.  A 
rectangular  hole,  circular  holes,  and  slats  were  used  to  set  as  a  vent  at 
the  RWE  end  plate.  Since  the  shape  and  number  of  holes  did  not  change  the 
resulting  waveform  at  Station  3,  the  station  of  most  interest,  it  is 
recommended  that  the  currently  used  end  flange  assembly  be  used  as  a  RWE 
directly  on  the  BRL  2.44  m  shock  tube.  Only  slight  changes  will  be  needed  for 
its  use. 

The  NASA-Ames  hydrocode  was  used  to  predict  several  shots  that  are  compared 
to  the  corresponding  experimental  shots.  The  predictions  compare  favorably  to 
the  experimental  data.  Pressure- time  records  are  shown  with  pertinent 
parameter  tables  for  use  of  a  rarefaction  wave  eliminator  at  the  BRL  2.44  m 
shock  tube. 

Appendix  B  shows  a  sketch  of  the  end  flange  assembly  that  might  easily 
be  adapted  for  use  as  a  rarefaction  wave  eliminator  by  using  either  a  solid 
plate  or  as  a  vented  plate.  The  small  l/48th  scale  shock  tube  was  found  to 
produce  waveforms  representative  of  the  large  2.44  m  BRL  shock  tube  at  all 
positions  where  pressure- time  records  were  obtained. 
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APPENDIX  A 


PRESSURE-TIME  RECORDS 
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Figure  A-l.  Pressure-time  records  for  solid  plate  RWE 
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Figure  A-2,  (Continued)  Pres sure- time  records  for  rectangular 
vented  RWE. 
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Figure  A- 2. 


(Continued)  Pressure-time  records  for  rectangular 
vented  RWE. 
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Figure  A-3.  Pressure-time  records  for  a  single  circular  hole 
vented  RWE 
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Figure  A-3.  (Continued)  Pressure-time  records  for  a  single 
circular  hole  vented  RWE 


94 


PRESSURE,  KPA 


PRESSURE,  KPA  PRESSURE,  KPA  PRESSURE,  KPA 


200 


TEST:  RUE,  2  IN  ST 
SHOT:  46 

STATION:  3  STAG  F 


TIME,  MSEC 


TEST:  RUE,  2  IN  ST 


TEST:  RUE,  2  IN  ST 


Figure  A-3.  (Continued)  Pressure-time  records  for  a  single 
circular  hole  vented  RWE 
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Figure  A-4.  Pressure-time  records  for  circular  vented  RWE  -  10  holes 


99 


TEST:  WE.  2  3N  ST 


TEST:  WE.  2  IN  ST 


TEST:  WE.  2  JN  ST 


TEST:  WE,  2  W  ST 


Figure  A-4.  (Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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Figure  A-4.  (Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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Figure  A-4.  (Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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Figure  A-4. 


(Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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Figure  A-4.  (Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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Figure  A-4.  (Continued)  Pressure-time  records  for  circular 
vented  RWE  -  10  holes 
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APPENDIX  B 

DRAWINGS  OF  ELIMINATOR  FOR  2.44  m  SHOCK  TUBE 
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Figure  B-l.  Side  View  RWE 


e  B-2.  End  View  RWE 
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Research  Inc. 

ATTN:  John  R.  Rempei 
2483  East  Bayshore  Road 
Palo  Alto,  CA  94303 

1  Effects  Technology,  Inc. 

ATTN :  E .  Anderson 
5383  Holister  Avenue 
Santa  Barbara,  CA  93105 

1  General  Electric  Co.  -  TEMPO 
ATTN:  DAS I AC 

816  State  Street;  Drawer  QQ 
Santa  Barbara,  CA  93102 

1  General  Electric  Co.  -  TEMPO 
7800  Marble  Avenue,  NE 
Suite  5 

Albuquerque,  NM  87110 

1  H-Tech  Laboratories,  Inc. 
ATTN :  B .  Hartenbaum 
P.O.  Box  1686 
Santa  Monica,  CA  90406 

1  Hughes  Aircraft  Company 
Systems  Development  Lab 
ATTN:  Dr.  A.  Puckett 
Centinela  8  Teale  Streets 
Culver  City,  CA  90230 
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Copies  Organization  Copies  Organization 


1  Ion  Physics  Corporation 
ATTN:  Technical  Library 
South  Bedford  Street 
Burlington,  MA  01803 

1  Kaman -Nuclear 

ATTN:  Dr.  D.  Sachs 

1500  Garden  of  the  Gods  Road 

Colorado  Springs,  CO  80907 

1  Kaman  Avidyne,  Division  of 
Kaman  Sciences 
ATTN:  Dr.  J.  Ray  Ruetenik 
83  2nd  Ave,  NW  Industrial  Park 
Burlington,  MA  01830 

1  KTECH  Corporation 

ATTN:  Dr.  Donald  V.  Keller 
911  Pennsylvania  NE 
Albuquerque,  NM  87110 

1  Lockheed  Missiles  8  Space  Co, 
Inc. 

Div  of  Lockheed  Aircraft  Corp 
ATTN:  J.  Nickell 
P.O.  Box  504 
Sunnyvale,  CA  94088 

1  Management  Science  Associates 
ATTN:  Kenneth  Kaplan 
P.O.  Box  239 
Los  Altos,  CA  94022 

1  Martin  Marietta  Aerospace 
Orlando  Division 
ATTN:  A.  Ossin 
P.O.  Box  5837 
Orlando,  FL  32805 

1  Maxwell  Laboratories,  Inc. 
ATTN:  A.  Kolb 
9244  Balbao  Avenue 
San  Diego,  CA  92123 


1  McDonnell  Douglas  Astronautics 

Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

2  Merrity  Cases,  Inc. 

ATTN:  J.  L.  Merritt 

Technical  Library 
P.O.  Box  1206 
Redlands,  CA  92373 

1  H.  L.  Murphy  Associates 
Box  1727 

San  Mateo,  CA  94401 

3  Physics  International  Company 
ATTN:  Document  Control 

Emil  Kovtun  (2  Cys) 

2700  Merced  Street 
San  Leandro,  CA  94577 

3  R8D  Associates 

ATTN :  Technical  Library 
Jerry  Carpenter 
Allen  Kuhl 
P.O.  Box  9695 
Marina  del  Rey,  CA  09291 

1  Sandia  Laboratories 
ATTN:  Dr.  J.  Kennedy 
Albuquerque,  NM  87115 

2  Science  Applications,  Inc. 

ATTN:  Joseph  McGahan 

Dr.  John  Cockayne 
1710  Goodridge  Dr.,  P.O.  Box  1303 
McLean,  VA  22102 
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No.  of 

Copies  Organization 

1  Systems,  Science  and  Software 
ATTN:  Technical  Library 
P.0.  Box  1620 
La  Jolla,  CA  92037 

1  Teledyne- Brown  Engineering 

Cummings  Research  Park 
Huntsville,  AL  35807 

1  Union  Carbide  Corporation 
Oak  Ridge  National  Lab 
ATTN:  Technical  Library 
P.O.  Box  X 

Oak  Ridge,  TN  37830 

1  Battelle  Memorial  Institute 
ATTN :  Technical  Library 
505  King  Avenue 
Columbus,  OH  43201 

1  Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  MD  20707 

1  Lovelace  Research  Institute 
ATTN:  Dr.  D.  Richmond 
P.O.  Box  5890 
Albuquerque,  NM  87115 

1  Massachusetts  Institute  of 
Technology 

Aerophysics  Laboratory 
77  Massachusetts  Avenue 
Cambridge,  MA  02139 

1  New  Mexico  Institute  of  Mining 
and  Technology 
ATTN:  Mr.  P.  McClain 
Socorro,  NM  87801 


No.  of 

Copies  Organization 

1  Northwestern  Michigan  College 
Traverse  City,  MI  49584 


1  Southwest  Research  Institute 
ATTN:  Dr.  W.  Baker 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1  Research  Institute  of  Temple 
University 

ATTN :  Technical  Library 
Philadelphia,  PA  19144 

1  Texas  Tech  University 

Dept  of  Civil  Engineering 
ATTN:  Mr.  Joseph  E.  Minor 
Lubbock,  TX  79409 

1  University  of  Arkansas 
Department  of  Physics 
ATTN:  Prof  0.  Zinke 
Fayetteville,  AR  72701 

1  University  of  California 
Lawrence  Livermore  Lab 
Technical  Library  Division 
ATTN:  Technical  Library 

Dr.  Donald  N.  Montan 
P.O.  Box  808 
Livermore,  CA  94550 

2  University  of  Denver 
Denver  Research  Institute 
ATTN:  Mr.  John  Wisotski 
P.O.  Box  10127 

Denver,  CO  90210 

1  J.  D.  Haltiwanger  Consulting Svcs 
B106a  Civil  Engineering  Bldg. 

208  N.  Romine  St. 

Urbana,  IL  61801 
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1  University  of  Maryland 
Department  of  Physics 
College  Park,  MD  20742 

1  The  University  of  New  Mexico 
Eric  H.  Wang  Civil  Eng'g  Res  Fac 
ATTN:  Technical  Library 
University  Station,  Box  188 
Albuquerque,  NM  87131 

1  University  of  Oklahoma 
Department  of  Physics 
ATTN:  Prof.  R.  Fowler 
440  W.  Brooks,  Rm  131 
Norman,  OK  73069 


Aberdeen  Proving  Ground 


Director,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Mr.  R.  Norman,  GWD 


Cdr,  USATECOM 

ATTN:  DRSTE-TO-F 


Dir,  USACSL 
Bldg  E3516,  EA 

ATTN:  DRDAR-CLB-PA 
DRDAR-CLN 
DRDAR-CLJ-L 
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USER  EVALUATION  OP  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


p 


Director 

US  Array  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BLA-S 

Aberdeen  Proving  Ground,  MD  21005 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE.  *300 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  12062  WASHINGTON.DC 
~ POSTAGE  WILL  BE  PAID  BY  DEPARTMENT  OF  THE  ARMY 


Director 

US  Army  Ballistic  Research  Laboratory 

ATTN :  DRDAR- BLA-S 

Aberdeen  Proving  Ground,  MD  21005 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


FOLD  HERL 


